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Abstract 

Mechanisms are elucidated underlying the existence of dynamical systems whose 
generic solutions approach asymptotically (at large time) isochronous evolutions: all 
their dependent variables tend asymptotically to functions periodic with the same 
fixed period. We focus on two such mechanisms, emphasizing their generality and 
illustrating each of them via a representative example. The first example belongs 
to a recently discovered class of integrable indeed solvable many-body problems. 
The second example consists of a broad class of (generally nonintegrable) models 
obtained by deforming appropriately the well-known (integrable and isochronous) 
many-body problem with inverse-cube two-body forces and a one-body linear ("har- 
monic oscillator") force. 
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1 Introduction 



Over the last three-four decades major progress occurred in the discovery and 
understanding of integrable dynamical systems with a finite or infinite num- 
ber of degrees of freedom, and over the last decade the possibility was noticed 
and exploited to identify and investigate many isochronous dynamical systems 
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characterized by a time evolution completely periodie (i. e., periodic in all de- 
grees of freedom) with the same period. This isochronous evolution might 
prevail in the entire (natural) phase space of the model under consideration 
(one talks then of an entirely isochronous system), implying of course that such 
a model is certainly integrable; or it might only prevail in an open (hence fully 
dimensional) region of its (natural) phase space, a phenomenology now known 
to characterize large families of nonintegrable dynamical systems possibly fea- 
turing quite complicated ("chaotic") behaviors outside the isochronous phase 
space region (for a review of these developments, see [1,2] ). In the present 
paper we discuss another, perhaps more interesting, phenomenology, namely 
dynamical systems whose generic solutions approach asymptotically (at large 
time) isochronous evolutions: all their dependent variables tend asymptotically 
to functions periodic with the same fixed period. The definition of such dy- 
namical systems is provided by the simultaneous validity of the two formulas 

lim [z„ (t) - 5„ (t)] = , n = l,...,N, (la) 

t — *-H-oo 

Zn{t + f) =Zn{t) , n= l,...,iV . (lb) 

Notation: the (generally complex^ but see below) numbers Zn (t) denote 
the N dependent variables of the dynamical system under consideration; we 
restrict consideration to the case when is a finite positive integer; the real 
variable t denotes the time; the functions Zn (t) characterize the asymptotic 
behavior of the dynamical system via (fTal) and the periodicity requirement 
(llb|) they satisfy characterizes the property of asymptotic isochronicity. This 
property is supposed to hold in an open (hence fully dimensional) region of the 
phase space of the dynamical system under consideration (possibly coinciding 
with its entire natural phase space): hence the dependent variables Zn (t) de- 
note here (the N components of) a generic solution of the dynamical system 
evolving (at least for sufficiently large time) within that region, while the func- 
tions Zn (t) , which shall generally be different for different solutions Zn (t) , are 
required to satisfy the periodicity property ffTbl) with the fixed period T (the 
same for all the solutions in the phase space region under consideration). Of 
course the formula (fTa|) does not define uniquely - for a given iV- vector z (t) 
- a corresponding A^- vector z_{t): the time-dependent A^- vector z_{t) is only 
identified by (ITa!) up to arbitrary corrections whose effects disappear in the 
asymptotic limit t oo. The property of asymptotic isochronicity is guar- 
anteed provided there exist just one A^-vector z(t) satisfying both relations 
([T]), for every generic solution z_{t) in an open, fully dimensional, region of 
phase space - namely for every solution z_{t) in that region of phase space, 
except possibly for some exceptional, generally singular, solutions belonging 
to a lower dimensional sector of that phase space region. 

The elementary idea underlying the identification of large classes of such 
asymptotically isochronous dynamical systems is to start from isochronous 
systems and then modify them by introducing a deformation whose effects 
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are significant through the time evolution yet disappear at large time: so that 
the modified systems loose their isochronous character (at finite times) but in 
some sense retain it (at large times) as the dominant feature characterizing 
their asymptotic behavior. 

There are several possible ways to implement this strategy in order to manu- 
facture asymptotically isochronous systems: some are rather trivial, some less 
so. This kind of judgement is of course subjective: for instance we tend to 
think that an important requirement for such systems to be deemed "interest- 
ing" is that they be autonomous - because the interest of dynamical systems is 
also related to their potential usefulness in order to model natural phenomena, 
which are generally described by autonomous evolution equations - and more- 
over because the freedom to introduce instead an explicit time dependence in 
the equations of motion of a dynamical system would provide too easy a way 
to infiuence more or less at will the asymptotic behavior of such a system. But 
of course the difference between autonomous and nonautonomous systems is 
unessential, since any nonautonomous system can be made autonomous by 
treating time itself as an additional dependent variable. 

In this paper we focus on two mechanisms yielding (autonomous) asymptot- 
ically isochronous systems, and illustrate each of them via a representative 
example. The first example (see Section 2) belongs to a recently discovered 
class of integrable indeed solvable many-body problems [3]; in this case we 
eventually focus on as simple and specific an example as possible, which is 
also suitable to exhibit some numerical results - but we trust our presentation 
is adequate to illustrate the generality of the approach. In this case the pe- 
riodic behavior prevailing asymptotically corresponds to a special solution of 
the dynamical system under consideration belonging to a region of phase space 
with positive codimension - albeit not an isolated solution of this system, so 
not quite identifiable as a limit cycle. Hence this model might be considered 
a representative example of a phenomenology characterized by the presence 
of some kind of friction. The second example (see Section 3) consists of a 
broad class of models obtained by deforming appropriately the well-known 
(see for instance [4]) integrable and isochronous one-dimensional many-body 
problem with inverse-cube two-body forces and a one-body linear ( "harmonic 
oscillator") force; the alert reader will again appreciate the generality of the 
approach, even though we illustrate it by focusing on a specific model (also 
restricting consideration to real dependent variables). In this second case the 
time-dependent A^-vector to which the solutions of the model tend asymptot- 
ically is not restricted to be in a sector of phase space with positive codimen- 
sion and is generally not itself a solution of the asymptotically isochronous 
A^-body model, so this phenomenology does not correspond to what is gener- 
ally referred to as a limit cycle behavior. In each of these two cases we back the 
qualitative understanding of the origin of the relevant phenomenology with a 
proof of its actual emergence, see ([T]). A section entitled "Outlook" in which 
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we elaborate tersely on the generality of this phenomenology concludes the 
paper. 



2 An asymptotically isochronous class of solvable many-body prob- 
lems 



A particular mechanism to manufacture integrable, indeed solvable, dynamical 
systems interpretable as many-body problems inasmuch as they are charac- 
terized by Newtonian equations of motion ("acceleration equal force") was 
introduced about three decades ago [5] and has been subsequently exploited 
to identify and investigate several such systems (for reviews of these devel- 
opments see for instance [2,4]). The idea is to exploit the nonlinear relation 
among the N coefficients Cm {t) of a (for definiteness, monic) time-dependent 
polynomial of degree N and its N zeros Zn (t): 



N N 



^ (z, t) = + ^ c„ it) =\[[z-z^ it)] . (2) 

m=l n=l 

A class of such systems is characterized by the fact that the N coefficients 
Cm {t) evolve in time according to a system of linear second-order constant- 
coefficient ODEs, the solution of which is a purely algebraic task (requiring 
essentially the diagonalization of an explicitly known matrix of order N) . The 
determination of the corresponding time evolution of the N zeros Zn {t) is 
therefore as well a purely algebraic task: computing the N zeros of a known 
polynomial. And it so happens that in many cases [2,4,5] this time evolution 
is indeed interpretable as that characterizing a Newtonian A-body problem - 
hence a solvable A-body problem, since its solution can be achieved by purely 
algebraic means. 

Indeed the solution z„ (t) of such a model is reduced to finding the A zeros of a 
polynomial of degree A in the {complex) variable z, see ([2]), whose coefficients 
Cm {t) generally evolve exponentially in time, typically 

N 

Cm (t) = E {T^^'^^^i^-^^ exp [A(^'+)t] + ^^^'-^ul^-^ exp [A^^'-)^] } , (3) 

where the 2A constants 7(^'='=) are arbitrary (to be determined by the initial 
data Zn (0) , i„ (0) in the context of the initial-value problem for the A-body 
system) and the 2A numbers A*^^'"*^-* respectively the quantities w^'^-* are the 
eigenvalues respectively the (components of the) eigenvectors of the matrix 
eigenvalue problem characterizing, as explained above, the dynamics of this 
system. Note that these eigenvalues and eigenvectors are associated to the 
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dynamical problem under consideration: they do not depend on the initial 
data identifying a particular solution, namely they are the same for all the 
solutions of the system. 

It is now clear (and indeed well known [2,4,5]) that if the 2N eigenvalues 
are all integer multiples of a single imaginary number iu (with > 0), 
_\(^>±) = ik^^'^'iuo with the 2N numbers k^^'^ arbitrary integers (positive or 
negative, but not vanishing), then the polynomial ijj{z,t) is clearly periodic 
with the (possibly nonprimitive) period 

T = ^, (4a) 

CO 

i^{z,t + T)^^{z,t) , (4b) 

hence all its zeros Zn (t) are as well periodic with this same period or possibly 
with a (generally small [6]) integer multiple p of this period, T = pT, due 
to the possibility that they exchange their role through the time evolution. 
Hence the corresponding A^-body problem is isochronous. 

And it is as well plain that if, out of the 2N eigenvalues A*^^'^\ only a (nonempty) 
subset have the property indicated above while all the others feature a neg- 
ative real part, then the many-body problem in question is asymptotically 
isochronous. This observation is not new, see for instance Section 4.2.3 of Ref. 
[4] (entitled "Some special cases: models with a hmit cycle, models with con- 
fined and periodic motions, Hamiltonian models, translation-invariant mod- 
els, models featuring equilibrium and spiraling configurations, models featur- 
ing only completely periodic motions"); but, to the best of our knowledge, 
this mechanism yielding asymptotically isochronous many-body problems was 
never analyzed in explicit detail (including the display of numerical results). 
This is what we do in this section, by focusing on a specific model whose inte- 
grable, indeed solvable, character has been ascertained only quite recently [3]. 



2.1 A specific example 



This A^-body problem (with iV > 3) is characterized by the Newtonian equa- 
tions of motion 



Zn — —a-iZn + 02-2, 



n 9 
4 



2a, 



zi + l 



zi 



2a4Zj, 



+2 



N 

E 

m=l,mjtn 



ZnZm + 02 + Qs^n + ^4 (4 " 1) 



n 



1,...,N 



(5a) 
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where the 4 "couphng constants" aj are a priori arbitrary complex numbers, 
superimposed dots denote time-differentiations and the rest of the notation is 
self-evident. The solvable character of this A^-body problem hinges [3] upon 
the following 4 restrictions on its initial data: 

E —A — = , E 2 = , (5b) 

which are then sufficient [3] to guarantee that, throughout the time evolution, 

implying that for this model it is justified to assume that only the evolution 
of — 2 particles is determined by the Newtonian equations of motion flSa|) . 
while the evolution of the remaining two is determined by these conditions, 
see (l5cll. 



Then the evolution of the N "particle coordinates" Zn [t) - taking generally 
place in the complex z-plane - coincides with the evolution of the zeros of 
a monic polynomial of degree in the variable z analogous to -0 {z, t) , see 
([2]), but more specifically reading as follows [3]: 

N-3 

Ip {Z, t) = TTn (z) + E [Cm (t) TTN-m (z)] + Cn (t) , (6a) 
m=l 

TTm{z) = z"' - e^—z^ - e^+imz, m = 0,1, N , (6b) 

£m = 1 if m is even , = if m is odd . (6c) 

And the coefficients Cm (t) evolve indeed according to formulas analogous to 
([3]), but more specifically reading as follows [3]: 



N 

(t) = {7^''+^«^^'+^ exp [A(^'+)t] + ^^''-^ui^-^ exp [X^'^-k] } , 

e=i,ejLN-i,N-2 

m = l,...,N — 3 and m = N , (7a) 



A(^'±) = ^i^, Aj = al + Ae[a, + i2N-e-3)a,] , 

e = l,...,N -3 ,N . (7b) 

Note that the coupling constant 03 does not appear explicitly in these formulas, 
but of course all 4 couphng constants Oj do play a role in determining the 
quantities m^'^^ appearing in the right-hand side of (17a|l . 
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We now restrict attention to the = 3 case, since this is sufficient, indeed con- 
venient, for exhibiting quite exphcitly an asymptotically isochronous model. 
Then the only relevant coefficient (see fl7al) ) is 

C3 (t) = 7+ exp (A+t) + 7_ exp (X^t) , (8a) 

A± = '""^^ ^ , A^ = al + 12a2 , (8b) 

where the somewhat simplified notation we are now using is we trust self- 
explanatory (and note that in this case with = 3 the eigenvalues X± only 
depend on the two coupling constants oi and 02). Correspondingly, the po- 
sitions of the 3 moving particles are the 3 zeros Zn (t) of the third-degree 
polynomial 

3 

i; {z, t) = TTs (z) + C3 (t) = z^-3z + C3 (t) = l[[z~Zn (t)] • (8c) 

n=l 

Note that these 3 zeros automatically satisfy the requirements f l5cl) . which 
corresponds [3] to the condition that the partial derivative of -0 {z, t) with 
respect to z vanish at z = ±1, ifj^ (±1, t) = 0. 

Assume now that the two coupling constants ai and 02 entail, via (l8b1) . 

X+ = iu , A_ = —a + i(3 , (9a) 

with a positive, a > 0, a; also positive, > (for definiteness) , and [3 real but 
otherwise arbitrary. This indeed happens provided 

uj {[3 + ia) , , 

ai = a-i{(3 + u) , 02 = ^^ -. (10) 

It is now plain that the asymptotic condition ( |TaI) holds now with Zn (t) being 
the three roots of the polynomial z^ — 3z + 7+ exp (iut) , 

3 

Z^-3Z + 7+ exp {iLUt) =l[[z-Zn (t)] , (11) 

n=l 

which provide of course also the special solution of the model (with A^ = 3) 
corresponding to initial data such that 7_ vanishes (see fl8al) ). And it is as well 
plain that the time evolution of this polynomial is periodic with period T, see 
(Hal) , hence the corresponding evolution of each of its 3 zeros is clearly periodic 
with periods T, 2T or 3T, depending whether that zero does not "exchange 
its role" through the motion with another zero or does so with one or with 
both the other two zeros. 

We complete this section by displaying one specific example, namely the so- 
lution of the system of ODEs with A^ = 3, a; = 27r implying T = 1 (see 
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Fig. 1. Trajectory of zi{t) in the complex z-plane from i = to t = 50 (see text) 

Re .zl 



fi I 




Fig. 2. Plot of Rezi as a function of t (see text) 

iajjj, as = a4 = 0, ai and 02 given by ffTOj) with a = 0.1 and /3 
with initial data 



-3, and 



Zl (0) = -2.1702823 + 0.180214311, ii (0) = 1.2487698 + 0.769412971, 
Z2 (0) = 0.71910399 - 0.891492881, Z2 (0) = -2.7507203 + 1.31025001, 
Z3 (0) = 1.4511783 + 0.711278571, Z3 (0) = 1.5019505 - 2.07966301, 



satisfying the conditions (l5b| and entailing 7+ = 0.5 + ?, 7_ = 3 — 3z (see (l8a|) ). 
The results displayed are, from t = to t = 50, the trajectory of Zi (t) in the 
complex z-plane (Fig. 1), the real part of zi (t) as a function of t (Fig. 2) (the 
behavior of the imaginary part is qualitatively analogous) and the evolution 
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D(t) 




40 t 



Fig. 3. Plot of the distance D{t) (see IH^ ) 

of the quantity 

D (t) = |C3 (t) - 7+ exp {%ujt)\ with C3 (t) = -z^ (t) (t) z^ (t) (12) 

(Fig. 3) that clearly provides a measure of the distance of this solution z (t) 
from its periodic limit z{t) (see ( IHcl) and (fTTl) . as well as ([I])). The numerical 
integration has been performed with an embedded Runge-Kutta method of 
order 8(5,3) with automatic step size control, as developed by Prince and 
Dormand [7]; the integration and the graphical output have been performed 
with the software Dynamics Solver developed by J. AguirregabiriaPI The 
results displayed have been obtained by integrating numerically the system of 
ODEs (jSD, checking throughout the integration the validity of the conditions 
(!5cl) as well as the two conditions 

Zl (t) + Z2 (t) + 23 (t) = , Z, (t) Z2 (t) + Z2 (t) Z3 (t) + Z; (t) Z^ (t) = -3 (13) 

(see (l8c|l ). The results reported are just a representative example of several 
numerical computations we did with different parameters and initial data, 
computations which were found to be quite reliable and stable unless the time 
evolution entailed a near collision of particles or their passage close to the 
special values z = ±1 (see fl5a|) ). 



This software is available at http://tp.lc.ehu.es/jma/ds/ds.htnil 
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3 A (generally nonintegrable) class of asymptotically isochronous 
many-body models 

In this section we consider a class of asymptotically isochronous models ob- 
tained by deforming the well-known integrable A'^-body problem with two-body 
inverse cube forces and a one-body linear force, which is of course isochronous 
when no deformation is present [4]. In particular we focus on the following 
equations of motion: 

1 ^ _ 

Xn + -u^Xn = g^ (Xn-Xm) ^ + F {w , X, x) , n = l,...,N , (14a) 

w = w[a\ogw — f {w,x,x)] , (14b) 

with 

< u;(0) < 1 . (14c) 

Here is an arbitrary positive integer (A^ > 2); the dependent variables 
Xfi = Xfi (t) may be interpreted as the coordinates of A^ particles evolving 
according to the Newtonian ("acceleration equal force") equations of motion 
f ll4al) ; these variables hereafter assumed to be all real (until we mention 

below to what extend the results change if the variables allowed to be 

complex), and x denotes of course the A^- vector with components Xn (this has 
motivated the notational replacement of the particle coordinates Zn with Xn, 
to be kept in mind when comparing the formulas written in this section with 
those written in the preceeding sections); likewise the auxiliary dependent 
variable w = w {t) evolves according to the first-order ODE fll4bp with an 
initial condition satisfying the inequalities fll4cp (but clearly, see below, one 
could replace this first-order ODE with an appropriate second-order "New- 
tonian" ODE); t denotes of course the (real) independent variable ("time": 
ranging from the initial time t = to the asymptotic time t = +00), and 
superimposed dots denote again differentiations with respect to this variable; 
u, and a are three positive (but otherwise arbitrary) constants; the main 
restriction on the, otherwise arbitrary, function F{w,x,v) is that it vanish 
when w vanishes, 

F{0,x,v)=0, (15a) 

and the main restrictions on the function / {w,x,v) is that it entail via ( ]14bl) 
a (very fast: see below) asymptotic vanishing (as t —>■ 00) of the auxiliary 
variable w (t) , 

lim \w{t)] = . (16a) 

t — >+oo 

A condition generally sufficient (but by no means necessary) to cause this is 
clearly (see fll4bl) with fll4cp and below) the requirement that f{w,x,v) be 
finite and nonnegative, 

0< fiw,x,v)<a\ (16b) 
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for all (real) values of w, x and v', it is indeed plain (for a proof, see below) 
that these conditions together with (]14b|) entail the inequalities 



0<w{t)<[w (0)]^"P("*) , 



(16c) 



hence (see fll4cl) and recall that a > 0) the auxiliary variable w (t) is always 
positive and vanishes asymptotically faster than exponentially, 



lim [w (t) exp (bt)] = 

t — >+oo 



(16d) 



with b any arbitrary constant. Restrictions on the dependence of the function 
F{w,x,v) upon the iV-vectors x and v are also required: a simple sufficient 
(but of course not necessary) condition, also encompassing fllSap . is that there 
exist a finite (positive) constant C and a positive number (3 such that 



\F {w,x,v)\ < C \w\ 



> 



(17) 



for all [real) values of w, x and v. Functions satisfying these conditions are 
for instance 



F {w,x, v) = Cw^ 
r {w,x,v) exp 



N 



n=l 



-1 



N 



n=l 



where An and Bn are arbitrary real constants. 

Our main result states that, for every (A^-vector) solution x (t) of this dynam- 
ical system, an (A^- vector) x (t) characterizing its asymptotic behavior (as 
t — i> +oo) via the formula (fTa|) (exists and) has the property to be completely 
periodic (i. e., periodic with the same period in all its components), see (llbp 
with T = T, see (l4al) . Of course this asymptotic A^- vector x{t) will depend 
on the solution x [t) under consideration - in particular, it will depend on the 
initial data, x (0) and x (0) , determining that solution in the context of the 
initial- value problem for the A^-body problem (|T^ : but let us re-emphasize 
that, for any arbitrary choice of these data (of course, satisfying the condition 
(0) 7^ Xm (0) for n ^ m, see f ll4al) ) it shall feature the property ([1]), namely 
all solutions x (t) of the system (HM shall feature the property of completely 
isochronous asymptotic periodicity ([T]) (with T = T, see f l4al) ) . 



This result is a natural consequence of the well-known fact (see for instance [4]) 
that all solutions of the system of Newtonian equations fll4al) without the F 
term in the right-hand side are completely periodic with period T, see (l4a|) . 
namely they all feature themselves the property (]lb|) with T = T. It stands 
therefore to reason that, if the function F {w,x,v) vanishes when w vanishes. 
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see fllSap . and if the time evolution (114bp of the auxihary variable w (t) entails 
that this dependent variable indeed vanishes asymptotically, see fll6al) . fast 
enough (see (116dp ). then asymptotically all solutions of our model (fn|) shall 
behave as the solutions of the same model without the F term, entailing the 
asymptotic phenomenology ([1]) with T = T, see fHaj) . 

To turn this hunch into a theorem a proof must be provided. This we do in the 
following subsection. Then in Section 4 we tersely discuss, again in the same 
qualitative vein as done above, to what extent the phenomenology described 
in this paper, and shown to occur in a specific, representative model, can be 
expected to occur in more general contexts. 



3.1 A theorem and its proof 

Theorem. The conditions (]16bl) and (fT7|) are sufficient to guarantee that every 
solution of the A^-body problem ffl^ with the three constants u, and a all 
positive yield the outcomes fll6al) and ([T]) with T = T, see ( Hal) : in particular 
they guarantee that there exists, corresponding to every solution x (t) of the 
A^-body problem f|T^ . an A^- vector x (t) satisfying both formulas ([1]) (of course, 
with Zn replaced by x„ and z„ by 

Proof. First of all let us prove the inequalities fll6cl) . obvious as they are. To 
this end we set 

(t) = [«; (0)pP[^W] , (18a) 

so that 

V? (0) = (18b) 

and (from fll4bl) ) 

ip{t)=a + f[w{t),xit),xit)]exp[-ip{t)] |log [w; (0)]r' , (18c) 

where we used the fact that log [w (0)] = — |log [w (0)]| , see (I14cl) . This ODE, 
together with the initial datum (]18b|) and the inequalities (116bp . clearly imply 
that (f [t) is positive and finite for < t < oo, indeed validity of the inequalities 

at < (p{t) < oo , < t < oo , (18d) 

which, via ffTSal) and ([He]), yield ffTGcj) . 

Next, let us introduce the counterpart of the Newtonian equations of motion 
f ll4al) . but without the F term in the right-hand side: 

Xn + -Uj'^ Xn = g"^ (Xn-Xm) ^ , n=l,...,N . (19) 

m=l,m^n 
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Here it is justified to use the notation x„ = x„ (t) for the dependent variables, 
since it is well-known [4] that all the solutions of this Newtonian A^-body 
problem are completely periodic with period T, see f l4al) . consistently with 
(fibl) with f = T. 



Let us now remark that, due to the strict positivity of g"^, this system of ODEs 
entails that 

\Xn (t) - Xm it)\ > C^ , > , H 171 , < t < OO , (20a) 

where (? is a time-independent constant that generally depends on the par- 
ticular solution under consideration but is certainly strictly positive, > 0. 
Likewise, again due to the strict positivity of g"^, the system of ODEs fll4al) 
with (^3) and (fTBcj) (entailing \F{w,x,v)\ < D, D = C \w {0)f) implies that 

|x„ (t) - (t)| > , > , nj^m, < t < oo , (20b) 

where is again a time-independent constant that generally depends on the 
particular solution under consideration but is certainly strictly positive, > 0. 
Moreover the systems of ODEs (1191) and (IT^ with (IT7|) and (I16cl) clearly imply 
that, for all (finite, positive) time, the functions x„ (t) and x„ (t) are finite. 

Let us now set 

e„ (t) = x„ (t) - 5„ (t) . (21) 

These functions C,n (t) satisfy - as implied by subtracting (fT9l) from fll4ap - 
the system of ODEs 

1 ^ 

+ + g'^ " ^rn] 'fnm {x, x) = F[w,X, x] (22a) 

m=l,m,y^n 

with 

H^nm (,21, XJ = - -3-— 1 -3 . (^//Dj 

Kp^n •^m) \-^n •^m) 

Note that the above bounds, (l20l) . as well as the finiteness of Xn and x„ for all 
(positive) time, guarantee that these functions Lpnm {x.x) remain finite for all 
time, namely that there always exist time-independent finite upper and lower 
bounds Lp± satisfied by them for all time, 

^~ < Vnm {X, X) < V9+ . (22c) 

These bounds depend of course on the particular solutions x and x under 
consideration, but let us re-emphasize that, for any such solutions, they are 
finite. 
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It is now clear that the theorem is proven if we can show that this system of 
ODEs admits a solution satisfying the asymptotic condition 



lim [^n{t)]=0, n=l,...,N (23) 

t — ^ + CX3 

(see f|Ta|) and fl2T]) ). As can be easily verified such a solution of ( !22|) is provided 
by the formula 

^nit)= dt'F[w{t'),xit'),xit')]Gn{t,t') , n=l,...,N, (24a) 



where the functions Gn {t, t') are the Green's functions associated with the left- 
hand side of the system of ODEs fl22al) . namely the solutions of the system of 
ODEs 



N 

+9^ E [Gn (t, t') - Cm (t, t')] ifnn, [x {t) , X {t)] = , t < t' 
m=l,m^n 

dGn (t, t' 



Gn {t, 0=0, 



dt 



= -1 , n = l,...,iV . (24b) 

t=t' 



Indeed, while these Green functions cannot be computed explicitly (since 
we do not know the A^- vectors x [t) and x {t) , hence neither the functions 
V^nm [2L{t) ,x{t)]), it is plain from the linear character of this system of ODEs 
and from the bounds fl22cl) that these Green functions can grow (in modulus) 
at most exponentially as t oo and/or t' —>■ oo ; so that the faster than ex- 
ponential vanishing of F [w {t') ,x{t') ,x(t')] as t' — > oo (implied by (ITTl) with 
(]16dp ) entails that the integral in the right-hand side of the solution formula 
fl24al) vanishes asymptotically (as t — > oo). □ 

Remark. It is clear how this example could have been made more general by 
allowing the function F appearing in the right hand side of fll4ap to depend 
on the index n, and/or by replacing the single auxiliary variable w {t) by a 
J- vector w (t) with J an arbitrary positive integer, and so on; without inval- 
idating our conclusion, but complicating our proof. Let us also re-emphasize 
that the hypotheses made above to prove this theorem are sufficient but by 
no means necessary for its validity. More specific, and possibly considerably 
less stringent, conditions yielding an analogous conclusion can and will be 
introduced whenever this kind of result shall be considered in specific (possi- 
bly applicative) contexts. Our motivation to assume here quite simple (hence 
overly stringent) hypotheses is because we are just interested to show that the 
main idea discussed in this paper does indeed work. □ 
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4 Outlook 



Clearly the kind of approaches illustrated above via the detailed treatment of 
two specific examples can be applied much more widely: it will be particularly 
interesting to do so in specific applicative contexts. 

A natural point of departure for such applications are isochronous systems, 
namely models whose generic solutions - in their entire natural phase space, 
or in open, hence fully dimensional, regions of it - are completely periodic 
(i. e., periodic in all their degrees of freedom) with the same fixed period 
(independent of the initial data, provided they stay within the isochronicity 
region). As recently pointed out (see for instance [2]), quite a lot of dynam- 
ical systems can be modified so that they become isochronous, entailing the 
conclusion that isochronous systems are not rare. Each of these isochronous 
systems can then be further extended - along the hues obviously suggested 
by the treatment detailed above, see in particular the specific case treated 
in Section 3 - in order to generate classes of asymptotically isochronous sys- 
tems, namely systems featuring open, hence fully dimensional, regions in their 
natural phase space (possibly including all of it) in which all (or almost all) 
their solutions display asymptotically a completely periodic behavior with the 
same fixed period, see ([1]). The technique to manufacture such generalized 
systems is clearly suggested by the examples treated above: of course these 
systems could be autonomous, as the examples treated above, or they might 
feature an explicit time-dependence, as could have been included in the sys- 
tem treated in Section 3 by assuming the functions F and / to also feature 
an explicit time dependence (but autonomous systems are generally more in- 
teresting than nonautonomous ones). 

Often the natural context to investigate isochronous systems is in the complex 
rather than the real [2,4] - although every system with complex dependent 
variables can of course be reformulated as a system with twice as many real 
dependent variables. Hence it may be of interest to mention how the findings 
detailed in Section 3 would be affected if the dependent variables Xn and w in 
the model (IT^ were allowed to be complex - keeping of course real the time t 
and positive the constant u, while the constant could now also be complex. 
It is then well known [2,4] that the isochronous character of the motions still 
prevails for the {integrable indeed solvable) many-body problem fll4al) with- 
out the F term (i. e., with an identically vanishing F; see (1191) ) - describing 
motions taking place in the complex 2;-plane rather than on the real line. But 
in the complex context the isochronous behavior is a bit different than in the 
real context: the phase space is then divided into sectors separated by lower- 
dimensional manifolds characterized by solutions which hit a singularity at a 
finite time due to a particle collision; an event forbidden in the real case with 
positive g^, when the particles move on the real axis and the two-body force. 
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singular at zero separation, is repulsive, see f ll4ap . but which can happen in 
the complex case, although not for generic initial data. In the different sectors 
the motion is still completely periodic, but with different periods, characteriz- 
ing each sector and being (generally rather small [6]) integer multiples of the 
basic period T, see fBal) . Accordingly, the generic solution of the (generally 
nonintegrahle) generalized model f|T^ will be nonsingular throughout its time 
evolution and it shall eventually settle within a sector, approaching asymptot- 
ically one of the completely periodic solutions in that sector of the {integrable) 
model (I14al) with identically vanishing F. 

A somewhat analogous outcome obtains for the model analogous to f|T^ but 
with f ll4ap replaced by 



featuring (A^ — 1) different coupling constants g^^^ acting among every par- 
ticle pair. In this case the model without F is generally not integrable, yet 
(if considered in the complex, namely without restricting the dependent vari- 
ables Zn - nor, for that matter, the coupling constants g^^ - to be real) it still 
does feature an open, hence fully dimensional, region in its phase space where 
all solutions are completely periodic with the same period T, see fHa|) [2,8]; 
while in other regions of its phase space it might also be periodic but with 
periods T = pT where the numbers p are integers but might be very large, or 
it might even display an aperiodic, quite complicated (in some sense chaotic) 
behavior [9] (for recent progress in the understanding of this phenomenology 
see [10-13]). It then stands to reason that the solutions of the generalized 
model (fT4l) with (]14ap replaced by ( !25l) (and of course x in ( ]14bl) replaced 
by z_) shall again approach asymptotically solutions - including, from open 
regions of initial data, completely periodic ones - of the model (1251) without F: 
entailing a remarkable, and quite rich, phenomenology. Clearly our motivation 
to mention this specific model is because of its prototypical role: indeed, the 
main aspects of this phenomenology shall also characterize the large class of 
isochronous (but by no means necessarily integrable) systems that can now 
be manufactured [2], once they are extended by adding to their equations of 
motion other, fairly general, terms having the property to disappear asymp- 
totically (as t +oo), as a consequence of the very dynamics implied by 
these extended equations of motion. 

In conclusion let us re-emphasize that these results (as indeed all mathemati- 
cally correct findings) might well be deemed remarkable or trivial, depending 
on the level of understanding of the reader. Once their foundation is under- 
stood, it becomes obvious how they can be extended to many other models - 
suggesting an ample applicative potential. But these developments exceed the 
scope of this paper. 
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